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ABSTRACT 

Near the minor axis of the Galactic bulge, at latitudes b < —5°, the red giant clump stars are 
split into two components along the line of sight. We investigate this split using the three fields 
from the ARGOS survey that lie on the minor axis at {I, b) = (0°, -5°), (0°, -7.5°), (0°, -10°). 
The separation is evident for stars with [Fc/H] > —0.5 in the two higher-latitude fields, but not 
in the field at b = —5°. Stars with [Fc/H] < —0.5 do not show the split. We compare the spatial 
distribution and kinematics of the clump stars with predictions from an evolutionary N-body 
model of a bulge that grew from a disk via bar-related instabilities. The density distribution of 
the peanut-shaped model is depressed near its minor axis. This produces a bimodal distribution 
of stars along the line of sight through the bulge near its minor axis, very much as seen in our 
observations. The observed and modelled kinematics of the two groups of stars are also similar. 
We conclude that the split red clump of the bulge is probably a generic feature of boxy/peanut 
bulges that grew from disks, and that the disk from which the bulge grew had relatively few stars 
with [Fe/H] < -0.5. 

Subject headings: Galaxy - abundances, kinematics, bulge, structure: Stars - late-type 



1. INTRODUCTION 

The Milky Way has a boxy/peanut-shaped 
bulge, as seen in the 2MASS star counts (Lopez- 
Corredoira et al. 2005) and the COBE/DIRBE 
near-infrared light distribution (Dwek et al. 1995). 
The long axis of this bulge/bar lies in the Galactic 
plane and points into the first Galactic quadrant 
at an angle of about 20° to the Sun-center line 
(Gerhard 2002). The distribution of stars as- 
sociated with the boxy /peanut structure shows 
complexity at higher Galactic latitudes. Near the 



minor axis, at b < —5°, recent analyses of op- 
tical and near-IR stellar photometry by Saito et 
al. (2011), Nataf et al. (2010) and McWilliam & 
Zoccali (2010) show that the red giant clump stars 
are split into two components separated by about 
0.65 mag along the line of sight. McWilliam & 
Zoccali (2010) argue that the split clump is as- 
sociated with an underlying X-shaped structure 
in the bulge. Such X-shaped structure has been 
noted earlier for other boxy/peanut bulges: (e.g. 
Whitmore & Bell (1988)). 
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We have recently completed the Abundances 
and Radial velocity Galactic Origins Survey (AR- 
GOS) which is a spectroscopic survey of 28,000 
stars, mostly clump giants, in 28 fields of the 
Galactic bulge and nearby disk. The ARGOS 
survey provides a clean sample of bulge giants 
with distances, radial velocities and metallicities. 
As part of this survey, we have data for about 
3000 stars in three two-degree fields along the 
minor axis, at {l,b) = (0°, -5°), (0°, -7.5°) and 
(0°,— 10°). These fields arc hereafter denoted as 
m0m5, m0m75 and mOmlO. Wc find an obvi- 
ous split in the magnitude distribution of the red 
clump stars along the minor axis in the higher- 
latitude fields, m0m75 and mOmlO. 

The objectives of this paper are (1) to measure 
the mctallicity range of the stars that show the 
split red clump near the minor axis of the bulge, 
(2) to determine whether the split red clump has 
an associated kinematic signature, and (3) to in- 
terpret these results in the context of an evolu- 
tionary N-body model. Our study complements 
previous photometric analyses of the stellar den- 
sity and magnitude distribution and expands upon 
previous spectroscopic studies: e.g Rangwala & 
Williams (2009); De Propis (2011). 

2. OBSERVATIONS 

The ARGOS observations were taken during 
42 nights in 2008-2011 with the Anglo Austrahan 
Telescope at Siding Spring Observatory. We used 
the AAOmega fibre-fed spectrograph (Sharp et al. 
2006) to observe about 340 stars simultaneously. 
The spectral region was 8400 — 8800 A at a resolu- 
tion of about 11,000. The locations of the 28 fields 
for the survey are shown in Figure [T] About 1000 
stars were observed in each field; the spectra have 
a typical signal to noise of 50 to 80 per resolution 
element. We selected stars from the Two Micron 
All Sky Survey (2MASS) in the magnitude range 
K = 11.5 to 14 with errors in J,K < 0.06 and 
with 2MASS flags set to reduce blends, contami- 
nation and stars of lower photometric quality. The 
colour-magnitude diagram (CMD) for our m0m75 
field in Figure [2] shows our colour and magnitude 
selection criteria. 

To make our selection cuts, we first dereddened 
the stars using the mean Schlegel et al. (1998) 
reddening for each field. For the three minor axis 
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Fig. 1. — The circles show the location of the 
ARGOS fields. The rectangle roughly indicates 
the outer limits of the Galactic bulge. The three 
minor axis fields are shown in red. The survey 
is now complete, except for the field at {I, b) = 
(20°,— 7.5°) for which only 600 stars were ob- 
served. 
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Fig. 2.— The K, (J - K) CMD for stars in the 
ARGOS m0m75 field at (/, b) = (0°, -7.5°). Stars 
for observation were selected from within the area 
outlined by the dashed box. which includes bulge 
giants with a wide range of mctallicity and ex- 
cludes most of the foreground disk dwarfs. The 
red clump giants lie near (J — A' )o = 0.65. 
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fields, the mean reddening is E{J—K) = 0.37, 0.15 
and 0.09 for our mOmS, m0ni75 and mOmlO fields 
respectively. Because the primary observations 
were made at the Ca triplet, and we wanted to 
observe stars of similar brightness together, we es- 
timated their /-magnitudes using a simple trans- 
formation from J—K to I—K derived from Bessell 
& Brett (1988). For the fibre allocations we then 
selected about 340 stars in each one-magnitude 
interval of /, from about 13 to 16. The integra- 
tion times for these three magnitude intervals were 
60 minutes, 120 minutes and 150 minutes, respec- 
tively. Full details of the ARGOS survey will be 
presented in Freeman et al. (2012, in preparation). 

3. DATA REDUCTION, RADIAL VE- 
LOCITIES AND STELLAR PARAM- 
ETERS 

We reduced our data using the standard 2dfdr 
pipelin(|3. This program subtracts bias and scat- 
tered light, extracts the individual stellar spectra, 
corrects for the flat field, calibrates wavelength 
and fibre throughput and subtracts sky. We se- 
lected the optimal sky subtraction option which 
gives much better sky subtraction, and used a 4th 
order polynomial for fitting the calibration lamp 
lines. 

We determined radial velocities and TV/, logf/, 
[Fe/H], and [a/Fe] for all stars, using our AR- 
GOS stellar pipeline (Ness et al. 2012, in prepara- 
tion), and estimated their distances via the Basti 
isochrones, assuming an age of 10 Gyrs. 

Radial velocities were obtained via cross cor- 
relation of the observed stellar spectra with syn- 
thetic spectra. The template spectra were gener- 
ated using the local thermodynamic equilibrium 
(LTE) stellar synthesis program MOOG (Sneden 
1973) using ID LTE atmospheres (Casteh & Ku- 
rucz 2003). 

The stellar parameters of the template spectra 
covered a range of [Fe/H] at a fixed gravity and 
temperature. The metallicity values used were 
[Fe/H] = -2.0, -1.0 and 0.0, for Te// = 5000 K 
and log g ~ 2.8 appropriate for a typical clump 
giant (Zhao et al. 2001). The best matching tem- 
plate to each observed star, as determined by the 
minimum error analysis via cross correlation, was 
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chosen to determine the radial velocity. The radial 
velocity errors at our observed resolution (11,000) 
and signal-to-noise 50 to 80 per resolution ele- 
ment) is < 1.2 km s^^. 

The heliocentric radial velocities were trans- 
formed to Galactocentric velocities. Wc took the 
circular velocity of the local standard of rest to be 
220 km s-i (Kerr & Lynden-BeU 1986) and the 
peculiar velocity of the sun to be 16.5 km s"^ in 
the direction {l,b) = (53°, 25°). 

Stellar parameters were determined using a 
minimum program written for this dataset; it 
compares each observed spectrum to a large grid 
of synthetic spectra. The grid of synthetic spec- 
tra was generated using the LTE stellar synthesis 
program MOOG 201(3 from the Castelh & Ku- 
rucz (2003) model atmospheres with no convective 
overshooting. The final grid covers a parameter 
range of T^// = 3750 - 6250 K, [Fe/H] = -5.0 to 
0.0, log .g = 0.0 to 5.0 and [a/Fe] = 0.0 to -hl.O. 
The grid steps are 250 K in T^ff, 0.5 in logg, 0.5 
in [Fe/H] (i.e. the steps of the original opacity 
distribution functions) and 0.1 dcx in [a/Fe], for 
a microturbulcnt velocity of 2.0 km s^^. 

The program was verified and tested using 
a number of calibration stars and fields and the 
Solar and Arcturus spectra from the Hinkle Atlas 
(Hinkle et al. 2000). FuU details of the pro- 
gram, including our linelist, grid, implementation 
and calibration will be given in a forthcoming pa- 
per by Ness et al. (2012, in preparation). 

Using our spectroscopic stellar parameters, we 
are able to select a clean sample of genuine clump 
giants. This sample can be further refined by using 
distances for our stars estimated from the derived 
stellar parameters. The apparent K magnitude of 
each star was corrected for interstellar extinction, 
using the Schlegel reddening at the position of each 
individual star and adopting Ak = 0.53i?(i? — V) 
(Schlegel et al. 1998). The correction is typically 
Ak — 0.25,0.1 and 0.05 for our three minor axis 
fields m0m5, m0m75 and mOmlO respectively. For 
the clump giants (identified as such from their 
temperatures and gravities) we adopted an ab- 
solute magnitude for the clump of Mk = —1.61 
(Alves 2000). For stars not in the clump, distances 
were estimated via fits to the Basti isochrones 
(Cassisi et al. 2006), adopting an age of 10 Gyr. 
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For the clump stars, the distance errors come from 
the uncertainties in (i) the interstellar reddening, 
(ii) the clump absolute magnitude and (iii) the 
2MASS apparent magnitude. We estimate our 
distance errors to be < 1.5 kpc. These distance 
estimates allow us to exclude stars which are not 
in the region of the bulge. 

In the following sections, we will discuss the 
magnitude distribution and kinematics of the split 
red clump seen in the bulge near its minor axis. 

4. THE ORIGIN OF THE SPLIT RED 
CLUMP 

McWiUiam & Zoccali (2010) found two red 
clump populations in the bulge at latitudes |6| > 
5.5°, separated by about 0.65 mag {Ko = 12.64 
and 13.29) or 2.3 kpc along the line of sight in 
their field near the minor axis {l,b) = (1°, —7.5°). 

From the three ARGOS survey fields on the mi- 
nor axis of the bulge, we have a sample of con- 
firmed red clump stars with known velocities, dis- 
tances and metallicities. First we would like to de- 
termine the metallicity range over which the split 
red clump is seen. The two higher-latitude fields 
m0m75 and mOmlO analysed separately show very 
similar results and to increase the number of stars 
in our analysis, they are combined for this discus- 
sion: the stars in these fields have b between —6.5 
and —11° and I between ±1°. To ensure a sample 
of clump stars lying in the bulge region, we select 
stars with log g = 1.9 to 3.1 and Ko = 12.38 to 
13.48. 

Examination of the K magnitude distribution 
of the clump stars function of [Fe/H] shows 
that the bimodality in the red clump magnitude 
distribution is very clear for the 271 red clump 
stars with [Fc/H] > —0.5 in the higher-latitude 
fields; it does not appear to be present for the 
206 more mctal-poor stars with —0.5 > [Fe/H] 
> —1.0. Figure |3] compares the magnitude dis- 
tributions for stars in the m0m75 -I- mOmlO fields 
across three [Fe/H] intervals: (i) [Fe/H] > 0, (u) 
> [Fe/H] > -0.5 and (in) -0.5 > [Fe/H] > -1.0. 
The binned data are shown in the upper panel 
and the normalised generalised distributions with 
Gaussian kernel a = 0.10 mag in the lower panel. 
In the generalised histograms, each star is repre- 
sented by a gaussian distribution with mean at 
the magnitude of the star and a = 0.10 mag (the 



maximum likely error of the K magnitudes for our 
stars) . We use generalised histograms to avoid the 
effects of arbitrary bin boundaries in conventional 
histograms. The generalised histograms clearly 
demonstrate the structure of the distribution. We 
do not show the stars with [Fe/H] < —1.0 here, 
because they probably belong to the inner metal- 
poor thick disk and halo (Ness et al., 2012, in 
preparation) . The peak of the magnitude distribu- 
tion for stars with —0.5 > [Fe/H] > — 1.0 lies at a 
fainter K magnitude than the mean of the more 
metal-rich stars, probably because we preferen- 
tially sample more of the more distant metal-poor 
stars in our survey cone at higher latitudes. The 
magnitude distribution of the stars with [Fe/H] 
< —1.0 is 15% narrower than for stars with —1.0 > 
[Fe/H] > —0.5 and, similarly, does not show the 
split clump. 

We note that in the higher-latitude fields the 
bimodality is seen very clearly for the most metal- 
rich red clump stars with [Fe/H] > 0, although 
they are kinematically different from the rest of 
the bulge. Their velocity dispersion is somewhat 
smaller (60 ± 5 km s~^) than for the more metal- 
poor bulge stars with [Fe/H] = -0.5 to 0.0 (81 ±3 
km s~^) which arc the dominant population in the 
bulge. The split for the stars with [Fe/H] > is 
marginally larger than for the stars with 0.0 > 
[Fe/H] > —0.5, and we sample more of the brighter 
metal rich stars on the near side of the bulge be- 
cause they lie preferentially at lower latitudes. We 
will return to this point in our discussion of the 
metallicity components of the bulge (Ness et al. 
2012, in preparation). 

We now examine the latitude dependence of the 
bimodality for the red clump stars which show the 
split red clump (those with [Fc/H] > —0.5). The 
left panels in Figure 0] compare the magnitude dis- 
tributions of the 377 stars in the higher-latitude 
combined (m0m75 + mOmlO) fields with the 290 
stars in the lower-latitude m0m5 field. The split 
red clump is very evident in the higher latitude 
sample, and appears to be absent in the m0m5 
sample. 

In order to interpret these observations, we 
compare them to an evolutionary N-body model 
of a boxy/peanut bulge from Athanassoula (2003) 
(A03). This model started as an isolated, axisym- 
metric galaxy with live disk and halo components, 
such that the disk is maximum (Sackett 1997; Ger- 
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Fig. 3. — Magnitude distributions for the clump 
stars in m0ni75 + mOmlO selected with 12.38 < 
iCo < 13.48, 1.9 < log g < 3.1 and (i) Fe/H > 
(81 stars) (ii) > Fc/H > -0.5 (190 stars), (iii) 
-1.0 > Fe/H > -0.5 (206 stars). The bimodal 
distribution is seen only for stars with [Fe/H] > 
-0.5 and is most prominent for stars [Fe/H] > 0. 
The top figure shows the binned data across the 
K magnitude range for each metallicity interval, 
and the bottom figure shows the normalised gen- 
eralised histograms (kernel a ~ 0.1) which show 
the minimum in the K magnitude distribution for 
stars [Fe/H] > -0.5. The vertical axis units for the 
bottom figure are arbitrary. 



hard 2006). The halo has a volume density de- 
scribed by: 
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where r is the radius, Mh is the halo mass, 7 
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Fig. 4. — The K-band magnitude distributions of 
clump stars red with [Fe/H] > -0.5 in the AR- 
GOS m0m5 field {upper left) and in the two com- 
bined fields m0m75 + mOmlO {lower left). Only 
stars with 12.38 < Ko < 13.48 are included. 
The magnitude distributions are shown as gener- 
ahscd histograms, using a Gaussian kernel with 
(7 = 0.10 mag. The panels on the right show 
the predicted magnitude distribution for the N- 
body model, adopting Mk = -1.51 for the clump 
stars. The vertical dashed lines show the magni- 
tudes of the two groups reported by Mc William & 
Zoccali (2010) in their bulge fields near the minor 
axis. The bimodal magnitude distribution is seen 
clearly for the bulge and the model in the higher 
latitude fields but only weakly in the m0m5 field. 
The vertical axis units are arbitrary. 

and Tc are the halo core and cut-off radii, re- 
spectively, and the constant a is given by a = 
[1 - 0r exp{q^) (1 - er/(5))]~\ where g = 7/rc 
(Hernquist 1993). 

The initial density distribution of the disc, for 
which there are 200 000 particles, is 



Pd{R, z) 



exp{-R/h) sech'^{ — ), (2) 



where R is the cylindrical radius, is the disc 
mass, h is the disc radial scale length and Zq is the 
disc vertical scale thickness. 
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We adopted zo/h = 0.2, r^/h = 10, Mh/Md = 5 
and 7/h = 5. The velocity dispersion in the disk 
was chosen such that Toomre's parameter is equal 
to 1.2 (see AOS). 

In this model, the bar which grows from the 
disk, buckles to form a boxy/peanut-bar/bulge 
(Athanassoula 2005). Although the model was not 
designed to match the Galactic bulge, it gives a 
good representation of its structure and kinemat- 
ics. 

The right panels in Figure |4] show the expected 
red clump magnitude distribution in the same 
fields, as derived from this model. The selection 
criteria for stars taken from the model were iden- 
tical to that of the observations. To select the 
roughly constant number of stars per unit magni- 
tude in the data, we scaled the model in (x, y, z) 
space by matching its projected exponential scale- 
height along the minor axis to the observed scale- 
height of the Galactic bulge. In velocity space, the 
model is scaled to match the velocity dispersion of 
the bulge stars in the ARGOS mOmS field. The 
model is then viewed from the Galactic plane at a 
distance of 8.0 kpc from the center and located so 
that the observer-center line is at 20° to the major 
axis of the bar. The absolute magnitude of the red 
clump stars in the model was set so that the mini- 
mum of the split red clump magnitude distribution 
in the (m0m75 -I- mOmlO) fields of the model is at 
Ko = 13.0, as observed for the ARGOS data in 
the same fields. The corresponding absolute mag- 
nitude of the red clump is Mk ~ —1.51. For com- 
parison, the Alves (2000) calibration is —1.65. For 
the model, the stars have been selected to mimic 
the roughly constant stars per magnitude in the 
data. 

The observed and model magnitude distribu- 
tions are qualitatively similar. The model mag- 
nitude distribution shows clear bimodality in the 
higher latitude fields (bottom right panel). In the 
lower mOmS field, the bimodality is almost absent 
(upper right panel). 

The orbit structure in the boxy/peanut bulge 
appears to generate a dip in the density distribu- 
tion near the minor axis at higher latitudes which 
produces the split structure in the red clump mag- 
nitude distribution. A related model by Shen et al. 
(2010) is similar. We have also examined several 
of Athanassoula's evolutionary N-body bar/bulge 
models, run with different initial conditions, and 



find that the split structure at higher latitudes 
near the minor axis appears to be a generic feature 
of bulges formed through disk instability. The de- 
tailed structure of the split (e.g. the magnitude 
separation of the two peaks) varies from model to 
model. 

We now visualize the split structure of the 
model in Figure [5l which shows the magnitude 
distributions in three slices of height above the 
Galactic plane. In the figure, (x, y) are coordi- 
nates projected on to the Galactic plane. The ori- 
gin is at the Galactic center, the y— axis lies along 
the Sun-center line, x increases towards decreas- 
ing longitude and z is the height above the plane. 
This figure is intended to illustrate the structure 
of the boxy/pcanut bulge. It is not intended for 
direct comparison with the data. It shows the den- 
sity distribution in slices of the z-height across the 
model parallel to the (x, y) plane. The three slices 
in z are 0.56 < z < 0.84, 0.91 < z < 1.20 and 
1.27 < z < 1.56 kpc. These slices correspond to 
the range of z for the three ARGOS fields on the 
minor axis (m0m5, m0m75 and mOmlO) where the 
line of sight intersects the y = plane. The figure 
shows only stars with \y\ < 3.5 kpc. The high- 
est slice 1.27 < z < 1.56 kpc shows a decrease 
in the density near the minor axis which we can 
associate with the boxy /peanut structure of the 
bulge. The intermediate slice 0.91 < z < 1.20 kpc 
also shows the two density maxima away from the 
minor axis. The lowest 0.56 < z < 0.84 kpc slice 
shows a central density peak with almost no sign 
of the double-peaked structure. This is consistent 
with the split clump results along the minor axis 
for the corresponding latitudes of the K magnitude 
distributions (Figure S]) where no split is seen at 
m0m5 and the split is seen at m0m75 -I- mOmlO,. 
The orbit structure of the boxy /peanut bulge ap- 
pears to generate a dip in the density distribution 
near the minor axis at higher latitudes, which pro- 
duces the split structure in the red clump magni- 
tude distribution. In projection, it would show the 
weak X-structure observed in other systems: e.g. 
IC 4767 (Whitmore & Bell 1988). 

In the Galactic bulge, it appears that this 
peanut structure involves only stars with [Fe/H] 
> -0.5. We conclude that the disk from which the 
bulge formed had relatively few stars with [Fe/H] 
< -0.5. This limit is unlikely to provide a useful 
guide on the time at which the bulge formation 
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Fig. 5. — Contour plots of the {x,y) density distribution of the N-body model for three slices in z: 0.56 < 
z < 0.84, 0.91 < z < 1.20 and 1.27 < z < 1.56 kpc corresponding to the three ARGOS bulge fields on the 
minor axis of the bulge. The sun lies at y = -8 kpc and the angle between the Sun-center line and the 
major axis of the bar is 20°. The higher slices show the dip near the minor axis in the density distribution 
associated with the structure of the bulge. The colorbar shows the normalised number of stars. This figure 
illustrates the peanut nature of the bulge model; it is not intended for direct comparison with the data. 



event occurred, because the chemical and dynam- 
ical evolution in the inner Galaxy were probably 
very rapid. 

The model seen in Figure [5] indicates that the 
split in the red clump should widen as the height 
above the plane increases. Figure [S] shows an in- 
creasing and more prominent separation between 
the two peaks in the density distribution with in- 
creasing height above the plane. For the stars with 
[Fe/H] > —0.5, which are involved in the split 
clump, we might therefore expect to see a differ- 
ence in the K magnitude distribution at m0m75 
compared to mOmlO. 

Figure [6] compares the K magnitude distribu- 
tion for the stars with [Fe/H] > —0.5 for the 
field m0m75 (143 stars) and the field mOmlO (147 
stars). The K-magnitude distribution for m0m75 
is clearly flatter than for mOmlO. At the higher 
latitude of mOmlO, the K magnitude of the fainter 
clump is fainter by about 0.06 mag leading to a 
slightly wider separation of the peaks for mOmlO. 
This small increase in separation corresponds to 
a linear increase in separation of 225 pc. When 
comparing Figures and [51 please note that the 
coordinates in Figure [5] are linear, while the hori- 
zontal axis of Figure |6] is in magnitudes. 

Figure[5]demonstrates the structure of boxy/peanut 
bulges in the central region. The flatter K- 




Fig. 6. — Magnitude distributions for the clump 
stars with [Fe/H] > -0.5 in m0m75 and mOmlO 
selected with 12.38 < Ko < 13.48, 1.9 < log g 
< 3.1. For m0m75 there are 143 stars and for 
mOmlO there are 147 stars. This figure shows that 
the clump separation is more clearly defined and 
is somewhat larger in K-magnitude for the higher 
latitude field. 



magnitude distribution for m0m75 and the deeper 
minimum and slightly wider separation of the 
peaks at mOmlO seen in our data appear consistent 
with the structure seen in the jzj- slices shown in 
Figure [S] For the data, taken in slices of latitude 
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rather than height below the plane, the distance 
from the plane increases with distance along the 
line of sight, and the relative density distributions 
of near and far clump stars along the line of sight 
are therefore affected. We investigate this below. 

The ARGOS data are not well suited to in- 
vestigating the relative numbers of brighter and 
fainter clump giants, because our stars were se- 
lected to be approximately uniform in their K- 
magnitude distribution. We do not use the AR- 
GOS data to test the decrease in density distri- 
bution of the total population of faint red clump 
compared to bright red clump along the line of 
sight function of latitude reported by Saito et 
al. (2011); McWiUiam & Zoccah (2010) and Nataf 
et al. (2010). We do see a bias in the apparent 
[Fe/H] distribution, with more metal-rich stars lo- 
cated in the near clump (see Figure 3) and more 
metal-poor stars at fainter magnitudes. This bias 
is associated with the metallicity gradient in the 
bulge (see Zoccah (2010), Ness & Freeman (2012)). 
The stars on the near side of the split clump are 
closer to the plane and are therefore more metal- 
rich in the mean, while the stars on the far side of 
the clump are at larger z— heights and are in the 
mean more metal-poor. 

The photometric studies by Saito et al. (2011), 
McWihiam & Zoccali (2010) and Nataf et al. 
(2010) examine the total density of the red clump 
stars in the central region. For the split red clump, 
these studies show that the near clump dominates 
and becomes brighter with respect to the faint 
clump at latitudes further from the plane. We 
can test for this effect in our model by exam- 
ining its stellar density distribution in slices of 
latitude rather than height. Figure [7] shows the 
density distribution of the stars in our model in 
the inner region in one-degree slices of latitude at 
|5| = 5°, 6°, 7° and 8°. Our model is symmetri- 
cal in z so we combine negative and positive lati- 
tudes. These one-degree slices are a subset of those 
mapped by Saito et al. (2011) at negative and 
positive latitudes, and are intended for compari- 
son with the Saito et al. (2011) maps. As in the 
Saito et al. (2011) data, the number of stars in 
the far clump in the models decreases relative to 
the near clump with increasing | latitude |, and the 
separation of the two clumps increases. Along a 
given line of sight, increasing distance corresponds 
to larger heights from the plane; at a latitude of 



|6| ~ 8°, the faint clump stars are about 300 pe 
further from the plane than the bright clump stars 
and contribute a significantly smaller fraction of 
the stellar density than the near red clump. This 
effect is a consequence of the decreasing density 
of each arm of the split clump as a function of 
distance from the plane. 

5. THE KINEMATICS OF THE SPLIT 
RED CLUMP 

We would like to know whether the stars in 
the near and far groups of the magnitude distribu- 
tion show different kinematic properties and, if so, 
whether these are reproduced by the model. First 
we examine the radial (line of sight) velocity dis- 
tributions of clump stars in the two latitude sam- 
ples, m0m5 and the combined m0m75 -(- mOmlO 
sample and compare them with the kinematics of 
the model. The minimum in the magnitude dis- 
tribution of the clump stars is seen at K = 13.0, 
so we divide the stars into a closer bright group 
(12.38 < K < 13.00) and a more distant faint 
group (13.00 < K < 13.48). 

The left panels of Figure |8] show the velocity 
distributions for the bright (solid curve) and faint 
(broken curve) groups in the two observed latitude 
samples. The velocity distributions are presented 
as generalised histograms, in which each stellar ve- 
locity is represented by a Gaussian with a = 35 km 
s~^ centered at its observed velocity; this avoids 
binning and suppresses high frequency noise in the 
distribution. This representation does not affect 
the mean velocities given below, which are cal- 
culated directly from the data: it is intended to 
make the form of the distribution more clear visu- 
ally. The numbers of stars in each subsample are 
shown in the panels. The right panels show the 
corresponding velocity distributions for the model. 
We recall that the model velocities were scaled so 
that the velocity dispersion of the model for the 
m0m5 sample agrees with the observed dispersion 
(97±4 km s~^). The velocity dispersion of the high 
latitude sample (77 ± 3 km s~^) is also similar to 
that of the model (70 ± 2 km s~^), with no signif- 
icant difference in dispersion between the bright 
and faint groups. For the low latitude field, the 
velocity offsets between bright and faint groups 
are —30 ± 12 km for the data and —24 ± 3 
km s~^ for the model. For the high latitude fields. 
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Fig. 7. — The (longitude, j/) density distribution for the N-body model in four one-degree-wide slices in 
latitude, at \b\ = 5°, 6° 7° and 8°. The sun lies at y = —8 kpc and the angle between the Sun-center line and 
the major axis of the bar is 20°. With increasing distance from the plane, the density of the more distant 
clump decreases relative to the nearer clump, and the separation between the nearer and more distant clumps 
increases. The colorbar shows the normalised number of stars. 



the offsets are +7 ± 9 km s^^ for the data and 
— 3±5 km for the model. The conclusion from 
this figure is that the velocity offsets between the 
near and far groups of the split clump are similar 
within the uncertainties for the data and model in 
both of the latitude regions. These velocity off- 
sets come mainly from streaming motions within 
the bar/bulge. The velocity offsets in the model 
are sensitive to the bar orientation relative to the 
Sun-center line; the velocity offsets agree with the 
observed offsets, within the errors, for a bar angle 
of 20°. 

In this section, we have considered only the 
kinematics of stars in ARGOS fields on the mi- 
nor axis. A detailed comparison of the kinematics 
for the model and the data for the other ARGOS 
fields, as shown in Figure 1, will be presented in a 
separate paper (Ness et al. 2012, in preparation). 

6. CONCLUSION 

The ARGOS survey provides a large sample of 
bulge stars with accurate abundances and radial 
velocities, relatively free of contamination from 
foreground and background stars. In this paper 
we used ARGOS stars near the minor axis of the 
bulge to investigate the properties of the split red 
clump. The split in the red clump magnitude dis- 
tribution is seen very clearly in the higher latitude 
fields {b = —6.5° to —11°) near the minor axis, for 



stars with [Fe/H] > —0.5. The more metal-poor 
red clump stars ([Fe/II]< —0.5) in the same re- 
gions of the bulge do not show the split, and it 
is almost absent at all metallicities for our m0m5 
field {b = -4° to -6°). 

We compared the structure and kinematics of 
the ARGOS stars with an N-body model of a 
bulge. In the model, the boxy/peanut bulge grew 
through the bar-forming and bar-buckling insta- 
bility of a disk of stars. Although the model was 
not constructed to fit the Galactic bulge, it is a 
good representation. The orbital structure of the 
boxy/peanut bulge leads to the central depres- 
sion in its peanut-shaped density distribution: see 
Athanassoula (2005). The central depression in 
turn leads to a split bimodal structure along the 
line of sight at higher latitudes which closely re- 
sembles the split observed in the ARGOS stars and 
in the Saito et al. (2011) data. At low latitude, 
the kinematics of the bulge stars in the bright and 
faint groups of the split red clump are also similar 
to those for the corresponding stars in the model. 

In projection, the density distribution would 
show an X-shaped structure, particularly after 
unsharp masking: see Whitmore & Bell (1988), 
Athanassoula (2005) and Bureau et al. (2006). It 
seems likely that the split red clump seen in the 
Galactic bulge is a generic feature of boxy/peanut 
bulges which grew from disks. We also conclude 
that the disk from which the bar/bulge formed 
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Fig. 8. — Velocity distributions of clump stars se- 
lected with Ko = 12.38 - 13.48. The central peak 
is stronger for the higher latitude and split clump, 
shown as generalised histograms using a Gaussian 
kernel with cr = 35 km s~^. Upper left: the ve- 
locity distribution for the m0m5 field. The dis- 
tributions for the near and far groups of the split 
clump are shown in the filled and dashed lines re- 
spectively. The numbers in the upper right corner 
of each panel show the number of stars in each 
group. Lower left: the velocity distribution for 
the two combined higher-latitude fields. The pan- 
els on the right show the velocity distributions of 
the same fields for the model. The normalised 
units on the vertical axis are arbitrary. 

would have had relatively few stars with [Fe/H] 
< -0.5. 
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A. APPENDIX 



An unexpected feature is seen in the velocity distribution of the bright group of clump stars on the minor 
axis of the bulge. Figure [S] compares the generalised histograms for the Galactocentric velocity distributions 
of stars in the bright group for the m0m5 field and for the higher latitude m0m75 and mOmlO sample. In 
this figure we use a much narrower kernel for the generalised histogram (ct = 2 km s~^. about double the 
measuring errors) to see the small scale structure. We can see this structure only because we have small 
measurement errors on our velocities of < 1.2 km s~^ which enables us to use this small kernel. Although 
the figure appears noisy because of the relatively small number of stars and the narrow kernel, both of 
the samples show a narrow and deep aligned minimum at a (heliocentric) radial velocity of 9 km s^^. No 
corresponding feature is seen in the model. 

To evaluate the probability of finding such an alignment of sharp deep minima by chance, we made a 
Monte Carlo simulation of pairs of samples. We are not testing whether the two distributions in Figure 
[5] are drawn fromthe same population: their different velocity dispersions show that they are not. We are 
estimating the probability that two random samples like those in the figure could have such deep aligned 
minima. In the simulation, the underlying distributions were gaussian with the same velocity dispersions as 
the two samples shown in Figure El We drew 1000 samples of 156 and 234 stars respectively from the two 
samples and found that aligned minima as deep as those seen in Figure |9] occurred with a probability of only 
0.007. Nevertheless this feature may not be real, and we draw attention to its presence without offering any 
suggestions for how it might come about. 
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Fig. 9. — Velocity distributions for stars in the bright group of the split clump, shown as generalised 
histograms using a Gaussian kernel with a — 2 km s~^. Upper panel: velocity distribution for the mOmS 
field. Lower panel: velocity distribution for the combined m0m75 and mOmlO fields. The numbers of stars 
in each subsample are shown in the upper right corner of each panel. The vertical dotted line at Vrad = 9 km 
s~^ shows the location of an apparent deep and narrow dip in the velocity distribution for both subsamples. 
No corresponding feature is seen for the stars in the faint group or in the model. The normalised units on 
the vertical axis are arbitrary. 
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